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Abstract 

Quantum walks have been used for developing quantum algorithms that 
outperform their classical analogues. In the context of spatial search algo- 
rithms, quantum hitting time plays an important role as the algorithm's 
stopping time. In this work, we analyze a decoherence model on Szegedy's 
quantum walk, based on percolation graphs. By performing averages over 
all possible evolution operators affected by the decoherence, we show that 
it is possible to define a decohcrent quantum hitting time and to prove that 
the quadratic speedup for the quantum hitting time and for the detection 
problem is still valid when the percolation probability is small. 

1 Introduction 

In Computer Science, Markov chains are employed in randomized algorithms 
such as searching algorithms on graphs. The expected time to reach a vertex for 
the first time, known as hitting time, plays an important role in those algorithms 
as the running time to find a solution. For instance, randomized algorithms are 
used to address the fc-SAT and the graph connectivity problem [T5] , 

Quantum Markov chains and quantum walks are the quantum analogue of 
classical Markov chains and random walks. There are at least three versions 
of quantum walks: 1) coined discrete-time quantum walks [I], 2) continuous- 
time quantum walks [8], 3) Szegedy's quantum walks [21]. All of them have 
been used for developing quantum algorithms that outperform their classical 
versions [191 IH [31 [TTJ [22] . These models use Hilbert spaces of different size 
and they seem to be non-equivalent. 

Using a new version of a discrete-time quantum walk model, Szegedy [21] 
showed that the quantum hitting time has a quadratic improvement over the 
classical one to detect a set of marked vertices for ergodic and symmetric Markov 
chains. This model was developed further by Magniez et al. [H], who described 
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a quantum algorithm with quadratic speedup for finding a marked vertex on re- 
versible, state-transitive Markov chains restricted to the case of only one marked 
vertex. Recently, Krovi et al. [11] showed that the quadratic speedup for finding 
a marked vertex also holds for any reversible Markov chain using a new interpo- 
lating algorithm. Santos and Portugal [TS] worked out the details of Szegedy's 
model on the complete graph, obtaining analytical results. 

When implementing quantum systems, decoherence problems are inevitable. 
It is impossible to completely isolate a real physical system and interactions with 
the environment will reduce or possibly destroy quantum coherence. These gen- 
erally undesired effects are present in quantum walk implementations. Hence, it 
is crucial to understand how decoherence affects them. Many papers address this 
issue. Brun et al. [5] analyzed the transition of a coined quantum walk to a clas- 
sical random walk by introducing decoherence on the coin operator. Kendon and 
Tregenna [10] also studied decoherence on the coin operator and highlighted some 
useful application of decoherence. Romanelli et al. [17] analyzed the decoherence 
produced by random broken links of a one-dimensional lattice. This technique 
was generalized for the two-dimensional case by Oliveira et al. [16]. Alagic and 
Russel [2] analyzed the effect of performing independent measurements on the 
continuous-time quantum walk on the hypercube. A review on decoherence on 
quantum walks can be found in [9]. 

Our contribution in this paper is to investigate a decoherence model on 
Szegedy's quantum walk inspired by percolation graphs. Our focus is on the 
quantum hitting time behavior in the presence of decoherence. By performing 
averages over all possible evolution operators affected by the decoherence, we are 
able to define a decoherent quantum hitting time. We show that the quadratic 
speedup on the quantum hitting time and for detecting a set of marked vertices 
for ergodic and symmetric Markov chains does not change significantly when the 
percolation probability is small. 

Decoherence inspired by percolation graphs was analyzed in many papers [17} 
[T6l 123 1 IT2 |, I13|. using the coined and the continuous-time quantum walk models. 
Previous to our work, decoherence on Szegedy's formalism was studied by Chiang 
and Gomez [6], who analyzed the sensibility to perturbation due to system's 
precision limitations, by adding a symmetric matrix E, representing the noise, 
to the transition probability matrix. The quadratic speedup vanishes when the 
magnitude of the noise \\E\\ is greater or equal to f2(<5(l — 6e)), where 5 is the 
spectral gap of the transition probability matrix of the graph and e is the ratio 
between the number of marked vertices and the number of vertices of the graph. 

This paper is organized as follows. In Sec. [21 we review Szegedy's quantum 
walk and the definition for the quantum hitting time. In Sec. [3] we describe 
the decoherence model on Szegedy's quantum walk and we define the decoherent 
quantum hitting time. The bound obtained for the decoherent quantum hitting 
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time is presented in Sec. HI In Sec. we analyze the Detection Problem and, in 
Sec. [6j we draw our conclusions. 



2 Szegedy's Quantum Walk 

Szegedy [21] has proposed a quantum walk driven by reflection operators in 
an associated bipartite graph obtained from the original one by a process of 
duplication. Let T(X,E) be a connected, undirected and non-bipartite graph, 
where X is the set of vertices and E is the set of edges. The stochastic matrix 
P associated with this graph is defined such that p xy is the inverse of the degree 
of the vertex x. Define a bipartite graph associated with T(X, E) through a 
process of duplication. X and Y are the sets of vertices of same cardinality of 
the bipartite graph. Each edge {xi,Xj} in E of the original graph F(X, E) is 
converted into two edges in the bipartite graph {xi,yj} and {yi,Xj}. 

To define a quantum walk in the bipartite graph, we associate with the graph 
a Hilbert space % n = H n <g> T~L n , where n = \X\ = \Y\. The computational basis 
of the first component is { cc) : x G X} and of the second : i/ £ Yj. The 
computational basis of is {la;, y) : x G X, y G Y}. The quantum walk on 
the bipartite graph is defined by the evolution operator Up given by 

Up:=n B K A , (1) 

where 

K A = 2AA T -I n 2, (2) 

Kb = 2BB T -I n 2, (3) 

2 2 

with the operators A : H n — > % n and B : % n — > H n defined as follows 

A = (4) 

yeY 

and 

|*x> = \x) ® I y\Jp^\v) I , (6) 




^y) = >. v^f) ®\y)- ( 7 ) 



In the bipartite graph, an application of Up corresponds to two quantum 
steps of the walk, from X to Y and from Y to X. We have to take the partial 
trace over the space associated with Y to obtain the state on the set X. 
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2.1 Coherent Quantum Hitting Time 



Instead of using the stochastic matrix P, Szegedy defined the quantum hitting 
time by using a modified evolution operator Up> associated with a modified 
stochastic matrix P', which is given by 

/ = / Pxy, X(£M- 

P *y \ 5 xy , x e M, W 

where M is the set of marked vertices. The initial condition of the quantum walk 
is 



W> = 7^v^M- (9) 

y&Y 

Note that |V>(0)) is an eigenvector of Up with eigenvalue 1. However, * s 
not an eigenvector of U pi in general. 



Definition 2.1 f21f The quantum hitting time Hp t M of a quantum walk with 
evolution operator Up given by Eq. |7p and initial condition k/>(0)) is defined as 
the least number of steps T such that 

F(T) > 1 — — , (10) 
n 

where m is the number of marked vertices, n is the number of vertices of the 
original graph and F(T) is 

T 

F(T) = > " |V>(0)}f, (11) 

t=o 

where Up, is the evolution operator after t steps using the modified stochastic 
matrix. 

The key operator to find the spectral decomposition of Upi is C = A T B |21| . 
The components of C xy are ^/p xy Pyx- We have to replace p xy by p' xy . Then, we 
have 

c -( F ; D- < i2 > 

where Pm is the matrix that we obtain from P by deleting its rows and columns 
indexed by M. The eigenvectors and eigenvalues of the evolution operator are 
obtained from the singular decomposition of the operator C. Therefore, its sin- 
gular values and vectors are directly related to the spectral decomposition of Pm 
and we can write each of these singular values as cos 9 for some < 9 < ^ . 
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Lemma 6 of Szegedy's paper [21] states that Hp t M is at most 

i-^v/T^aT <13) 



n 



where A' l5 . . . , X' n _ m are the eigenvalues of Pm and \v[), . . . , \ v' n _ m ) are the as- 
sociated normalized eigenvectors. Pm is the matrix obtained from P by leaving 
out all rows and columns indexed by some x € M. Coefficients are defined 

such that |u) = Ylk=T v k\ v 'k) f° r p) = — 7=1 j where 1 is the (n — m)-dimensional 

vector with entries equal to 1. 



3 Decoherent Quantum Hitting Time 

Percolation model was introduced in the context of a porous media for analyzing 
liquid flows, and is now a paradigm model of statistical physics [20J. The de- 
coherence model inspired by percolation can be explained in the following way: 
Suppose that a walker is on a vertex of a graph. Before moving to the neigh- 
boring vertices, each edge can be removed and each non-edge can be inserted 
with probability p. With probability 1 — p, an edge stays unchanged and the 
same is valid for a non-edge. After this change in the graph topology, the walker 
moves following the dynamics of the model. The original graph is reset and the 
process is repeated over in the next step. This decoherence model is called bond 
percolation. 

The occurrence probability of a given Pj is determined as follows. If < p < 1, 
then Pr(Pi) = (1 — p^ ac ~ ad p ad , where a c = n ^ n ~ 1 ^> [ s the number of edges of the 
complete graph with n vertices and is the number of edges removed plus the 
number of edges included to obtain Pj from P. If p = 0, Pr(Pj = P) = 1, and 
Pr(p / P) = 0. And, if p = 1, we have Pr(P = P) = 1, and Pr(P / P) = 0, 
where P is the complement of P. The evolution under maximum decoherence 
occurs when p = 1/2, since at each step a random graph is selected. 

Another decoherence model can be analyzed at this point. If only removal of 
edges are allowed (no insertions), a c must be replaced by the number of edges of 
the original graph and ad will be the number of removed edges. 

The dynamics under decoherence has a new behavior, because at each step 
the graph changes. This process changes the transition probability matrix, which 
also changes the evolution operator. Therefore, instead of having an usual walk 
evolving as \ip(t)) = Up\ip(0)), we have 

\m) = UkUp^ ■■■C/ Pl |V(0)) =: C/pJV(0)), (14) 

where p = (Pi, . . . , Pt-i, Pt) and Up t = Up t Up t _ 1 • • •Up 1 . Pi's are not necessarily 
equal and they are independent and identically distributed. They are obtained 
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from P and for each Pj we have a P/ associated depending on the cardinality of 
M. 

In this context, it is useful to define an operator that will represent the be- 
havior of the operators affected by the decoherence. Let 

U dec :=Y,Pr(P)Up, (15) 
P 

be the operator obtained by performing an average over all possible evolution 
operators affected by the decoherence. The following result shows that the av- 
erage over all possible sequences P, with size T, according to its probability 
distribution, is equal to Uj ec - 

Lemma 3.1 For t < T we have 

52Pr(Pr)U Pt = U t dec . (16) 

Pt 

Proof Since Pr(P T ) = ]jj =1 Pr(P), we have 

T 

2 Pr(PT)Up t Up t _ L ---Up, = PHPiWpUp^ ---Up, 

P T Pt i=1 

= E E • • E flI Pr (^)) U P* U r*-i -V* f ^Pr(Pi)U Pl J 

P T Pt-i P2 \i=2 / \ Pi J 

= E E • • • E Pr(Pr)Pr(Pr-i) ■ ■ ■ Pr(P t+1 )U t dec 

P T Pt-i Pt+i 
— u dec 



To define the quantum hitting time for the decoherent evolution, we have to 
do an average over all possible sequences P. Define, 



F dec (T) :=£Pr(P T ) ( ^— £ 



Lemma 3.2 

T 

2 



\up t \n > 2 1 



F dec (T) = 2 - — — £ <^(0)|^ ec |^(0))- (18) 

i=0 
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Proof Expanding Eq. (|18p and using that the initial condition and the evolution 
operators are real, we obtain 

F dec {T) = £Pr(P T ) (^^( 2 -2^(0)|^J^(0)>) ) 



Using Lemma l3.1( we obtain Eq. ()18|) . 

Now, we can naturally define the decoherent quantum hitting time (DQHT), 
using the expression of F dec obtained in Lemma 13.21 

Definition 3.3 The decoherent quantum hitting time Hp^M of a quantum walk 
with evolution operator Up given by Eq. f7J) and initial condition IV'(O)) is defined 
as the least number of steps T such that 

F dec (T)>l--. (20) 
n 

Note that when p = 0, we have the original definition, since U dec = Up/. 



4 Bounds on DQHT 

It is important to mention that we are considering ergodic Markov chains with 
symmetric transition matrix. Thus, the following result generalizes the result 
from Szegedy [21] by introducing a decoherence term. 

Theorem 4.1 The decoherent quantum hitting time Hp e \i of a quantum walk 
with evolution operator Up, given by Eq. (CP, initial condition \ip(0)), and p < 



l 

300a c £ 



where 



E i _ m £ arccos(Ap ' (21) 



is at most 



k> 



n—m o „.„ /n—m o 

v% 942 a c p sr^ uf. 



m / j 



k-- 
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Proof Using expression ([22]) and replacing \' k by cos 6^ , define 



n—m 2 I n—m 2 \ 

T = 8Y- V , k + 942a cP y V , k . (23) 

^ (l-e)Vl-cos^ \ti (l-e)Vl-cos^y 1 ; 

Using 1 — cos a > 2a 2 /5 (a £ (0, 7r/2]), we obtain 

n—m 2 

T - 13 ^ (Tr^r + 2355a ^( y.Tr^r) • (24) 
Using 1 — cos a < a 2 , we obtain 

n—m 2 /n—m 2 

T>8^ 7T ^ + M2a cP [y 7T ^) . (25) 




- (l-e)9 k 

Using £ = Ylk=i (i-e)e k ' we obtain 

8£ + 942a c p£ 2 < T < 13E + 2355a c ;p£ 2 . (26) 
Let us write the initial condition as |V>(0)) = |^a/-l) + IV'm)) where 



\1>m±) = s/v^\x)\y), (27) 

\4>m) = -7= Yl VPw\ x )\y)- ( 28 ) 



n 

yGJC 

Using Eq. ([!]), we obtain that ^m 1 -) = Ylk=T u kA\vk), where \v k ) are the eigen- 
vectors \v' k ) increased with m zeros to have the correct dimension, compatible 
with the dimension of C. Part of the singular decomposition of C is associated 
to the eigenspectrum of Pm- When Pm is symmetric, the singular vectors \w k ) 
are equal to ±k>fc), its sign depends on the sign of the associated eigenvalue of 
Pm- The vector |V>m x ) belongs to the two-dimensional space spanned by A|u>fc) 
and P>\v k ) with associated singular value cos(#fc). In this space, the evolution 
operator Upi is a sequence of two reflections around different axes. Then, we 
obtain a rotation with an angle that is twice of the one between the axes. 

We want to show that F<iec(T) is greater than or equal to 1 — ^ when T is 
in the range (f26j) . Using Eq. (fl~9j) and |V>(0)) = \iPm- 1 -) + h/'iw), we can write 
Fdec{T) as 

F dec (T) = 2 - 2{G M + G MM ± + G M ±), (29) 
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where 



1 T 

g m = E Pr ( p T) E (M u p t \M, (30) 



1 T 

Cm,mx = ^yE^^E (W^J^m) + (V'mI^J^)) (31) 



T 

G M ± 



1 T 

— E P < P t) E \Up t \ip M x). (32) 



T+ _ 



Let us establish bounds for Gm, G m m j_ and G^: 

T 



1 T 

G M < y— E Pr (^) E (V>m|Vm) = e. (33) 



Pt t=o 



Expression (131|) can be expanded in two terms 

T 

-Pr(Pv = (P 1 P'\\ ' 

T + 



1 T 

P T ^(P',...,P') *=° 

(34) 

The first term of Eq. (j34"|) is zero because k/'Af- 1 ) is in the space spanned by ^4|u>fc) 
and .B|ffc), which are eigenvectors of Up>. Taking e < 1/2, 

(^A/^l^pJ^A/) + (^M\Up t \ip M ±) < 2max{<V'M|V ; A/}, (^Af-L \^M±)} 

= 2 max{e, 1 — e} 

= 2(1-6). 

Then, we have 

G MjM x < 2(1 - e)(l - (1 - p) acT ) < 2(1 - e)a cP T. (35) 
Finally for G M ± , we have 

1 T 
G M ± = ¥ —Pr(PT = (P',...,P'))^2^MA U P'\^M^ + 

t=o 



1 T 



T+ _ 

P T ^(P',...,P') *=o 
/.. w-. -, a y n—m 2 T 

< ( f~ p)c £ E cos ») + a - - a -p) acT y 

k=l 6 i=0 

From Eq. (13) of Szegedy's paper [21], we know that 

1 n—m 2 T ^ n—m . 2 

1 y y cos{2te ) <7 J—Y 7 -i\ r . (36) 

k=l t=0 k=l y ' K 

Then, 

G A/ x < (1 - e) ( (1 - a cP T) ^ + o&t) + 0(p 2 ). (37) 



From Eqs. (J33J), ([35J), and d37J), we obtain 

Cm + C M)M ± + G A/ x < e + (1 - e) f ^ - 4a c pE + 3o c pT^ + 0(p 2 ). (38) 
Using that T is in the range (|26j) and p < 30Q * E , we have 

Gm + G MM ± + G M ± < — y~ . (39) 

Then, 

F dec (T) > 2 - 2 (^] = 1 - e. (40) 



Corollary 4.2 T/ie decoherent quantum hitting time Hp e %i of Up with respect 
to any M C X with m < n/2 and < p < 30Q * E , where 



n—m o 

E 



— Y 



^^arcooB(Ai)' 

is in O == = V where X(Pm) is the largest eigenvalue of Pm in absolute 

value. 

Proof Using 1 — cos a we obtain 



1 Z / 1 

E >ly-^-J — 1 — (4i) 

- 2 1-eV 1-cos^ 
fe=i 
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and 



P < 



(42) 



By substituting Eq. (|42|) for the expression of Hp e %i given by (|22|) and since 




(43) 



we conclude that H^ r is in O , 

Expression (|22p of Theorem 14.11 shows that the decoherent hitting time has 
an additional term that is proportional to the square of the usual term. If p is 
small enough, the contribution of the new term for the hitting time scales as a 
linear function in terms of p. Corollary 14.21 describes a range of p such that the 
quadratic speedup is valid. 



5 The Detection Problem 

Let M C 2 X be a set of non-empty subsets of M. The Detection Problem stands 
as the problem of finding out if the set of marked vertices is either empty or 
belongs to M.. 

Theorem 5.1 Assume that T is an upper bound for 

n-\M\ 2 /n-\M\ 2 \ 2 

"SttV^HStM' (44) 

where M runs through all elements of Ai (\' k and v k depend on M) and p obeys 
the same inequality of Theorem \4-l\ Then, the Detection Problem can be solved 
within time T with bounded two-sided error. 

Select < t < T uniformly random and let Ui, ■ ■ ■ ,Ut be unitary operators 
of the dynamics under decoherence. Algorithm [T] creates the state 

§|0> (|^(0)) + UtUt^.-.U^))) + i|l> (|V'(0)> - UtUt-x.-Uxl^O))) , (45) 

which has an additional control register of one qubit. Note that, when using 
Algorithm [H we do not know if Up or Up> are being used. By measuring the 
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Algorithm 1: Detect if marked 



input : Ui, ■ ■ ■ ,U t 

output: (M = 0) or 1 (M G M) 

begin 

prepare the state: 0) ^(0)); 
apply H to the control register; 

apply C(Ui), ■ ■ ■ , C(Ut), where C(U) is the controlled-U; 
apply H to the control register; 
measure in computational basis; 
if the control register is 1 then 
return 1 

else 

return 



state given by Eq. (|45p in the computational basis, the probabilities of having 
the control register in state and 1 are 



p(0) 

p(i) 



|V(0)) + £Wt_i...£/i|V(0)) 

\1>(o)) -u t u t -i...uMo)) 



(46) 
(47) 



Now, analyzing the algorithm after measuring (Step[6j), if M = 0, the control 
register is in state |o) with probability at least (1 — p) acT , that is the probability 
of having U = Up. If M G Ai, the control register is in state 1 1) with probability 
at least 



TT^TTT ^ Pr{Pr) ^ 1 1 1 VJ(0) > " ^ ^ (0) 
1 ' s t=o 



2 1 , 



(48) 



Eq. (|48p is obtained from Eq. (j47j) by performing two averages: an average on 
time, because the algorithm chooses a time t at random, and another average on 
the possible sequences Pt since the algorithm can be affected by the decoherence. 
Therefore, using that ^ < ^, we obtain 1 in the control register with probability 
at least |, which means that there is at least one marked element. 

This result can be improved if we consider the decoherence model which allows 
only removal of edges from the graph (insertions are not allowed). We can solve 
the Detection Problem within time T with bounded one-sided error, because the 
initial condition will be invariant under the action of any Up i when M is empty. 
In this case, the probability to obtain when there are no marked vertices is 1. 
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6 Conclusions 



We have proposed a decoherence model on Szegedy's quantum walk inspired by 
percolation graphs. This model is characterized by the possibility of removing 
or inserting edges at each time step with probability p. The graph probability 
matrix and the evolution operator change at each time step. We were able to 
define a decoherent quantum hitting time by using a new operator, which is 
obtained by performing an average over all possible evolution operators affected 
by the decoherence. Note that when the percolation probability p is zero, the 
evolution operator is equal to Szegedy's original definition and the decoherent 
hitting time is equal to the original quantum hitting time |21j . 

We have proved that, for p small enough, the decoherent quantum hitting 
time has an additional term which depends linearly on p, preserving the quadratic 
speedup on the quantum hitting time over the classical case for some range of 
p. Moreover, the Detection Problem can be solved with time of the order of the 
DQHT with bounded two-sided error, and with bounded one-sided error if the 
decoherence model does not allow edge insertions. 

It is interesting to analyze the behavior of the hitting time in the complement 
of the original graph. This case occurs when p = 1 and Udec = Upi. Since the 
initial condition is associated to the original graph P, it is easy to show that 
|^(0)) is invariant under the action of Upi. This occurs because the initial 
condition is a superposition over all edges of the original graph and those edges 
do not exist in its complement. Thus, (V'(O)) is an eigenvector of Upi with 
eigenvalue 1 and the quantum hitting time goes to infinity when p goes to 1 . 

In future works, it is interesting to obtain a better upper bound for the DQHT. 
The approximations used in Eqs. ([30]) to (f32|) help the analytical calculations, 
but seems to weaken the final value of the bound. It is also interesting to analyze 
the decoherence effects on spatial search algorithms. 
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